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Abstract 

We present a measurement of CP-violation parameters in the KgKgKg decay based 

on a sample of 275 x 10® BB pairs collected at the T(4S') resonance with the Belle detector 
at the KEKB energy-asymmetric e"*“e“ collider. One neutral B meson is fully reconstructed in 
the decay B^ and the flavor of the accompanying B meson is identified from its 

decay products. CP-violation parameters are obtained from the asymmetry in the distributions 
of the proper-time interval between the two B decays: S = -1-1.26 ± 0.68(stat) ± 0.20(syst) and 
A = -rO.54 ± 0.34(stat) ± 0.09(syst). 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 


*on leave from Nova Gorica Polytechnic, Nova Gorica 
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In the Standard Model (SM), CP violation arises from the Kobayashi-Maskawa phase ^ 
in the weak-interaction quark-mixing matrix. In particular, the SM predicts to a good 
approximation that S = —^/sin20i and M = 0 for both b ccs and b sqq transitions, 
where S (M) is a parameter for mixing-induced (direct) CP violation [3, ^/ = 1) 

corresponds to CP-even (-odd) hnal states, and (f)i is one of angles of the Unitarity Triangle. 
Measurements of time-dependent CP asymmetries in P° —>■ J/'ll ) O and related dec^ 
modes, which are governed by the b —> ccs transition, by the Belle and BaBar p 

collaborations already determine sin 20i rather precisely; the present world average value is 
sin201 = +0.726 ±0.037 0. 

CP-violation parameters in the flavor-changing 6 —> s transition are sensitive to phe¬ 
nomena at a very high-energy scale 0, 0. Belle measurements Q in the decay modes 
P° —>• (j)Kg, 0P°, K~^K~K^, fo{980)Kg, rj'Kg, ujKg, and which are dominated by 

the b —>■ sqq transition, yield sin20i = PO.dS^oji when all the modes are combined. Mea¬ 
surements by the BaBar collaboration also yield a similar deviation [l, 1I|- To elucidate the 
difference in CP asymmetries between b —> sqq and b —>• ccs transitions, it is essential to 
examine additional modes that may be sensitive to the same 6 —»• s amplitude. 

The P° decay to KgKgK^, which is a = ±1 state, is one of the most promising modes 
for this purpose l^l- Since there is no u quark in the hnal state, the decay is dominated by 
the b —> sss transition. Its branching fraction P(P° —>• K^K^K^) = (4.2l};3 ± 0.8) x 10“® 
was reported by Belle Q. In this Letter, we describe a measurement of CP asymmetries 
in the P° —> KgKgKg decay. 

In the decay chain T(4S') —P°P° —> /cp/tag, where one of the B mesons decays at 
time tcp to a CP eigenstate fcp and the other decays at time ttag to a hnal state /tag that 
distinguishes between P° and P°, the decay rate has a time dependence given by 


g-+4/TsO 

P(At) = —-- \ 1 ± g ■ S sin(AmrfAt) 


Atbo 


+A cos(AmrfAt) 


( 1 ) 


Here Tpo is the P° lifetime, Arud is the mass diherence between the two P° mass eigenstates. 
At is the time diherence tcp — ttag, and the 6 -havor charge is g = ±1 (—1) when the tagging 
B meson is a P° (- 8 °). 

At the KEKB energy-asymmetric e’''e“ (3.5 on 8.0 GeV) collider Q , the T (45*) resonance 
is produced with a Lorentz boost of 07 = 0.425 along the 2 ; direction which is antiparallel to 
the positron beamline. Since the P° and P° mesons are approximately at rest in the T (45*) 
center-of-mass system (cms). At can be determined from the displacement in z between the 
fcp and /tag decay vertos: At ~ {zcp - +ag)/(07c) = Az/(07c). 

The Belle detector [13 is a large-solid-angle magnetic spectrometer that consists of a 
silicon vertex detector (SVD), a 50-layer central drift chamber, an array of aerogel thresh¬ 
old Cherenkov counters, a barrel-like arrangement of time-of-hight scintillation counters, 
an electromagnetic calorimeter (ECL) and an iron hux-return instrumented to detect 
mesons and to identify muons. A 2.0 cm radius beampipe and a 3-layer SVD (SVD-I) were 
used for a 140 fb“^ data sample, while a 1.5 cm radius beampipe, a 4-layer silicon detector 
(SVD-II) and a small-cell inner drift chamber were used for an additional 113 fb“^ data 
sample. In total 275 x 10® BB pairs were accumulated. 

We reconstruct the P® —>• KgKgKg decay in the Kg~Kg~Kg~ or Kg~Kg~Kg^ hnal 
state, where the 7 r’'' 7 r“ (tt^tt®) state from a Kg decay is denoted as Kg~ {Kg^). Pairs of 
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oppositely charged tracks with the tt+tt” invariant mass within 0.012 GeV/c^ (~ 3 ct) of the 
nominal mass are used to reconstruct Kg~ candidates. The tt+tt” vertex is required to 
be displaced from the interaction point (IP) by a minimum transverse distance of 0.22 cm 
for candidates with p > 1.5 GeV/c and 0.08 cm for those with p < 1.5 GeV/c, where p 
is the momentum of K^. The angle in the transverse plane between the momentum 
vector and the direction dehned by the vertex and the IP should be less than 0.03 rad 
(0.1 rad) for the high (low) momentum candidates. The mismatch in the z direction at the 
vertex point for the two charged pion tracks should be less than 2.4 cm (1.8 cm) for the 
high (low) momentum candidates. After two good Kg~ candidates have been found which 
satisfy the criteria given above, looser requirements are applied for the third Kg~ candidate. 
The requirement on the transverse direction matching is relaxed to 0.2 rad (0.4 rad for low 
momentum candidates), and the mismatch of the two charged pions in the z direction is 
required to be less than 5 cm (1 cm if both pions have hits in the SVD). 

To select candidates, we reconstruct vr® candidates from pairs of photons with 
> 0.05 GeV, where is the photon energy measured with the EGL. The reconstructed 
7r° candidate is required to have an invariant mass between 0.08 and 0.15 GeV/c^ and mo¬ 
mentum above 0.1 GeV/c. candidates are required to have an invariant mass between 
0.47 and 0.52 GeV/c^, and a £t is performed with constraints on the vertex and 
masses to improve the invariant mass resolution. The candidate is combined with 
two good Kg~ candidates to reconstruct a meson. 

To identify —>• KgKgKg decays, we use the energy difference AE = 
and the beam-energy constrained mass Mbc = (A'b™')^)^ — (p™®)^, where is the beam 

energy in the cms, and and are the cms energy and momentum of the reconstructed 
B candidate, respectively. The B^ meson signal region is dehned as |Ai7| < 0.10 GeV for 
B^ Kj-Rj-Rj-, -0.15 < AE < 0.10 GeV for 5° ^ Rj-Rj-Rf, and 5.27 < 
Mbc < 5-29 GeV/c^ for both decays. To suppress the qq continuum background 

{q = u, d, s, c), we form a signal over background likelihood ratio TZg/h by combi ning 
likelihoods for two quantities; a Fisher discriminant of modihed Fox-Wolfram moments (l7[, 
and the cosine of the cms B^ hight direction.The requirement for TZg/h depends both on the 
decay mode and on the havor-tagging quality; after applying all other cuts, this rejects 94% 
of the qq background while retaining 75% of the signal. 

If both B^ —>■ Rg~Rg~Rg~ and Rg~Rg~R^ candidates are found in the same event, we 
choose the B^ —> Rg~Rg~Rg~ candidate. When multiple B^ Rg~Rg~Rg~ candidates 
are found, we prioritize those with three good Rs~ candidates. If more than one candidate 
still remain, we select the one with the smallest value for X](AM^+-)^, where AM^+- is 

the difference between the reconstructed and nominal mass of Rs~- For multiple B^ 
Rg~Rg~R^ candidates, we select the Rg~Rg~ pair that has the smallest X](AM^+-)^ 
value and the R^ candidate with the minimum of the constrained £t. 

We reject RgRgRg candidates if they are consistent with B^ —>• XcoKg {RgRg)Rg or 
B^ —>■ D^Rg —> {RgRg)Rg decays, i.e. if one of the Rg pairs has an invariant mass within 
±2ct of the Xco mass or mass, where a is the RgRg mass resolution. 

Figure^shows the Mbc and AE distributions for the reconstructed B^ —>■ RgRgRg candi¬ 
dates after havor tagging. The signal yield is determined from an unbinned two-dimensional 
maximum-likelihood £t to the AF-Mbc distribution. The Rg~Rg~Rg~ signal distribution 
is modeled with a Gaussian function (a sum of two Gaussian functions) for Mbc (AF). For 
B^ —> Rg~Rg~Rg^ decay, the signal is modeled with a two-dimensional smoothed his¬ 
togram obtained from Monte Garlo (MG) events. For the continuum background, we use 
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FIG. 1: (a) Mbc distribution within the AE signal region, (b) AE distribution within the Mbc 
signal region. Solid curves show the fit to signal plus background distributions, and dashed curves 
show the background contributions. 




18| for Mbc and a linear function for AE. The fits after flavor 
Kt~Kt~ events and 16±8 5° ^ Kt~Ki 


~Kg^ events for 


the ARGUS parameterization 
tagging yield 72 ± 10 

a total of 88 ± 13 ^ KgKgKg events in the signal region, where the errors are statisti¬ 

cal only. The obtained purity is 0.56 for the Kg~Kg~Kg~ and 0.40 for the Kg~Kg~Kg^ 
channels. We use events outside the signal region as well as a large MC sample to study 
the background components. The dominant background is from continuum. The contam¬ 
ination of —> XcoKg events in the —>■ sample is small (less than 2.6% at 

90% C.L.). The contributions from other BB events are negligibly small. The influence of 
these backgrounds is treated as a source of systematic uncertainty in the CP asymmetry 
measurement. Backgrounds from the decay B^ —> are found to be negligible. 

The 6-flavor of the accompanying B meson is identified from inclusive properties of par¬ 
ticles that are not associated with the reconstructed B^ —>■ K^KgKg candidates. The 
algorithm for flavor tagging is described in detail elsewhere US We use two parameters, 
q dehned in Eq. and r, to represent the tagging information. The parameter r is an 
event-by-event, MC-determined flavor-tagging dilution factor that ranges from r = 0 for no 
flavor discrimination to r = 1 for unambiguous flavor assignment. It is used only to sort 
data into six r intervals. The wrong tag fractions w for each of the r intervals and their 
differences Aw for B^ and B^ decays are determined from data 3 . 

The vertex position for B^ —> KgKgKg decays is obtained using Kg~ trajectories and 
a constraint on the IP; the IP prohle {ax — 100 fim, cXy ~ 5 /im, ~ 3 mm) is convolved 
with hnite B^ flight length in the plane perpendicular to the 2 ; axis. To reconstruct the Kg~ 
trajectory with sufficient resolution, both charged pions from the Kg decay are required to 
have enough SVD hits; at least one layer with hits on both sides and at least one additional z 
hit in other layers for SVD-1, and at least two layers with hits on both sides for SVD-11. The 
reconstruction efficiency depends both on the Kg~ momentum and on the SVD geometry. 
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The vertex efficiencies with SVD-II (86% for Kg~Kg~Kg~ and 74% for Kg~Kg~Kg^) are 
higher than those with SVD-I (79% for Kg~Kg~Kg~ and 62% for Kg~Kg~Kg^) because 
of the larger outer radius and the additional detector layer. The typical vertex resolution is 
about 97 /iui (113 /ini) for SVD-I (SVD-II) when two or three Kg~ candidates can be used. 
The resolution is worse when only one Kg~ can be used; the typical value is 152 pm (168 pm) 
for SVD-I (SVD-II), which is comparable to the /tag vertex resolution. The determination 
of the vertex of the jVgffiial state is the same as the analysis, and is described 

in detail elsewhere |l0|,l20|. 

We determine S and A by performing an unbinned maximum-likelihood £t to the observed 
At distribution. The probability density function (PDF) expected for the signal distribution, 
'Psig(At; S, A, q, w, Aw), is given by Eq. ((T)) after incorporating the effect of incorrect flavor 
assignment. The distribution is convolved with the proper-time interval resolution function, 
Rsig, which is a function of event-by-event vertex errors. 

We hnd from MC simulation that universal i?sig parameters used for measurements of 
CP asymmetries in the J/%IjK^ and related decays approximately describe the 

resolution for the 5° —>■ KgKgKg decay. To account for differences between Kg~ trajectories 
and charged tracks, additional parameters that rescale vertex errors are introduced. When 
only one Kg~ is used in the vertex £t, these parameters are determined from a fit to the At 
distribution of B^ J/'ll) candidates, where only the and the IP constraint are used 
for the vertex reconstruction. The procedure is the same as that for the B^ decay 

and is described elsewhere jl^- For events with two or three Kg~ used in the vertexing, 
we also hnd from MC simulation that the resolution is well described by the same i?sig 
parameterization with an additional correction function that depends on the number of 
Kg~ decays used for the vertex reconstruction. The form of this correction function is 
determined from a study using MC simulation. 

We determine the following likelihood value for each event v. 


Pi = (1 - foi) j /sigPsig(AF)^sig(Ati - At') 

+ (1 - /sig)Pbkg(AfOi?bkg(Ati - At') d{At') 
+ fo\Po\{Ati), 


( 2 ) 


where Poi is a broad Gaussian function that represents an outlier component with a small 
fraction /oi The signal probability /sig is calculated on an event-by-event basis from the 
function obtained by the AP-Mbc two-dimensional £t used to extract the signal yield. A 
PDF for background events, Pbkg, is modeled as a sum of exponential and prompt compo¬ 
nents, and is convolved with a sum of two Gaussians Pbkg- All parameters in Pbkg and Pbkg 
are determined by the £t to the At distribution of a background-enhanced control sample, 
i.e. events outside of the AP-Mbc signal region. We £x r^o and Am^ at their world-average 


values j21|. In order to reduce the statistical error on A, we include events without ver¬ 
tex information. The likelihood value in this case is obtained from the function of Eq. m 
integrated over At*. 

The only free parameters in the final fit are S and A, which are determined by maximiz¬ 
ing the likelihood function L = Yii PilAtf, S , A) where the product is over all events. An 
unbinned maximum likelihood £t to the 167 B^ —> KgKgKg candidates, containing 88 ± 13 
KgKgKg signal events, yields 


S = -1-1.26 ± 0.68(stat) ± 0.20(syst), 
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FIG. 2; Raw asymmetry in each At bin with (a) 0 < r < 0.5 and (b) 0.5 < r < 1.0. The 
solid curves show the result of the unbinned maximum-likelihood fit. The dashed curves show 
the SM expectation with (5, .4) = (—0.73,0). The numbers of candidate events used for the S 
measurement are 65 for (a) 0 < r < 0.5 and 52 for (b) 0.5 < r < 1.0. 


A = +0.54±0.34(stat) ±0.09(syst). 

We define the raw asymmetry in each At interval by {N+ — N^)/[N^ + N^), where fV+(_) is 
the nnmber of observed candidates with q = 1). The raw asymmetries in two regions 

of the flavor-tagging parameter r are shown in Fig.|21 Note that these are simple projections 
onto the At axis and do not reflect other event-by-event information (snch as the signal 
fraction, the wrong tag fraction and the vertex resolntion), which is in fact used in the 
unbinned maximum-likelihood £t for S and A. 

The systematic error is primarily due to the resolution function (±0.12 for S and ±0.04 
for A), the background fractions (±0.10 for S and ±0.03 for A), £t bias (±0.08 for S and 
±0.05 for .4.), and background modeling (±0.08 for S and ±0.01 for A). Other sources of 
systematic error are uncertainties in the wrong tag fraction (±0.04 for S and ±0.01 for A), 
physics parameters Am^ and r^o (±0.01 for S and ±0.01 for A), the vertexing (±0.02 for 
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S and ±0.05 for A), and the effect of tag side interference [23 (±0.02 for S and ±0.02 for 
A). We add each contribution in quadrature to obtain the total systematic errors. 

Various cross-checks of the measurement are performed. We reconstruct 
decays without using the charged kaon for the vertex reconstruction and apply the same £t 
procedure. We obtain = ±0.33lo!45(stat) and AkOko^k± = —0.32 ± 0.28(stat), 

which are consistent with no CP asymmetry. MC pseudo-experiments are generated to 
perform ensemble tests. We hnd that the statistical errors obtained in our measurement are 
all consistent with the expectations from the ensemble tests. We apply the same procedure 
to the P° —>■ J/tjjKg sample without J/'ijj daughter tracks for vertex reconstruction. We 
obtain = ±0.68 ± O.lO(stat) and Aj/^^o = ±0.02 ± 0.04(stat), which are in good 

agreement with the world average values 0]. We conclude that the vertex resolution for the 
P° —> KgKgKg decay is well understood. 

We use a frequent ist approach 0 to determine the statistical signihcance of the deviation 
from the SM. From 1-dimensional conhdence intervals for S, the case with S = —0.73 for 
P° —>• KgKgKg is rulod out at a 99.7% conhdence level, equivalent to 2.9a signihcance for 
Gaussian errors. 

In summary, we have performed the measurement of CP-violation parameters in the 
P° —>■ KgKgKg decay based on a sample of 275 x 10® BB pairs. The decay is dominated by 
the b ^ s havor-changing neutral current and the KgKgKg hnal state is a CP eigenstate. 
It is therefore sensitive to a possible new CP-violating phase beyond the SM. The result 
dihers from the SM expectation by 2.9 standard deviations. 

We thank the KEKB group for the excellent operation of the accelerator, the KEK 
cryogenics group for the efficient operation of the solenoid, and the KEK computer group 
and the Nil for valuable computing and Super-SINET network support. We acknowledge 
support from MEXT and JSPS (Japan); ARC and BEST (Australia); NSEC (contract 
No. 10175071, China); DST (India); the BK21 program of MOEHRD and the CHEP SRC 
program of KOSEF (Korea); KBN (contract No. 2P03B 01324, Poland); MIST (Russia); 
MHEST (Slovenia); SNSF (Switzerland); NSC and MCE (Taiwan); and DOE (USA). 

Note added .—As we were preparing to submit this paper, we became aware of a paper 
from the BaBar collaboration [2^ which reports on the branching fraction and CP asym¬ 
metries in the R® —>■ K^KgKg decay. 
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